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\^ ' Abstract 

We compute the leading contribution to the effective Hamihonian of 
SU{N) matrix theory in the hmit of large separation. We work with a 
gauge fixed Hamiltonian and use generalized Born-Oppenheimer approx- 



QQ ' imation, extending the recent work of Halpern and Schwartz for SU{2) 

I 1^. The answer turns out to be a free Hamiltonian for the coordinates 

along the flat directions of the potential. Applications to finding ground 
state candidates and calculation of the correction (surface) term to Witten 
index are discussed. 
Oh! 

' 1 Introduction 



We study the maximal supersymmetric gauge quantum mechanics. This model 
^ was first considered in papers Q, ||^, Later this Hamiltonian found appli- 

ed ' cations in the physics of supermembranes Q , j|] and _D-particles j6| . There is a 

remarkable conjecture, called Matrix theory, that the same Hamiltonian gives a 
nonperturbative description of M-theory |Q. This conjecture implies the exis- 
tence of a unique bound state at threshold for each SU{N) gauge group where 
N > 2 (though the existence of the state was conjectured earlier). 

The SU{2) case was recently studied in a number of papers ||^, |p^ , 
[ p] . In it was proved, using the computation of Witten index, that at least 
one ground state exists in SU (2) model. Some progress also have been achieved 
in the general SU{N) case | |T^ . In jT^ the principal contribution to 

Witten index was computed (see also ]14| and ]15||). Calculation of the surface 
contribution to the index was performed in jl^ under certain assumptions. The 
present paper justifies these assumptions. In |^ an asymptotic form of a wave 
function was proposed as a candidate ground state. 
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The present paper extends the ideas of Q and ||Tl| to higher A^'s. We 
consider the model in the asymptotic region where the coordinates along the 
flat directions of the potential become large. Our method is generalized Born- 
Oppenheimer approximation introduced in The authors of |^ tackle the 
problem of gauge invariance by working with a complete set of gauge invariant 
variables. It would be rather hard to implement this idea for higher iV's. In- 
stead, we apply the gauge fixing procedure of The result of our computation 
is that in the leading order the dynamics of flat coordinates is described by the 
free Hamiltonian. Possible potential term vanishes due to cancellations between 
many different terms. This result is not surprising, it is generally expected in 
supersymmetric theories. It is in agreement with SU{2) calculations of |pl| , 
. In the last section of the paper we discuss possible applications of the main 
result. 



2 The model and preliminaries 

In this section we introduce the model and do some preliminary technical work 
needed for the Born-Oppenheimer approximation. Namely, we write down the 
(partially) gauge fixed Hamiltonian, split it into 4 basic parts: free Hamiltonian 
for slow (Cartan) degrees of freedom, bosonic and fermionic oscillators, and the 
interaction part. After that we explain the quantization of the oscillators part. 
We refer the reader to the seminal paper Q for all the details of the gauge fixing 
procedure. 

The Hamiltonian of matrix theory is that of the lOD super Yang-Mills theory 
dimensionally reduced to 0+1 dimensions (in Aq = gauge). One can write it 
in the following form 

11 i 

^ = + ■^fABC<l>%<l>cfADE4>D(l>E " ■^fABC'I'A^Ba^ap^CP (1) 

where 0^, tt^ are real bosonic variables and Aao are real fermionic variables. 
The lower capital Latin indices are that of the adjoint representation of a real 
compact Lie algebra g with totally antisymmetric structure constants /abc- 
Denote by G the Lie group of g. The indices fi, i' — 1, ... ,9, a = 1, ... ,16 cor- 
respond to vector and real spinor representations of the group spin(9) respec- 
tively. The S0{9) gamma matrices F'^ are assumed to be real and symmetric. 
The canonical commutation relations between the variables are 

4>b] = -iS^^^AB , {^Aa, Ab/j} = 5ap5AB- 

From now on we set G — SU{N) that corresponds to the Matrix theory with the 
"center of mass" degrees of freedom being excluded. The Lie algebra g = su{N) 
consists of the traceless antihermitean N x N matrices (f). An invariant, positive 
definite inner product on g is defined as (0i,(/'2) = ^2rr((/)i(/>2). The Cartan 
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subalgebra consists then of diagonal matrices (f) G g. Let matrices Ta, A = 
1, . . . , — 1 be an orthonormal basis in su{N) with respect to this inner product 
such that the indices A = 1, . . . , N — 1 correspond to the Cartan subalgebra. 
Then fABC = {.TA,[TB,Tc]). 

Now we will briefly explain the gauge fixing procedure. Given any element 
4> G g there exist a unique matrix D such that (f> = UDU^^ for some U E G 
and D is of the form 



D = 









A2 



V 

where A„ are real numbers such that 



\ 





Xn J 



(2) 



N 

71 = 1 



0, 



Ai > A2 > . . . > A 



N- 



(3) 



The transformation U is defined up to multiplications by an arbitrary element of 
the Cartan subgroup. Thus, it is clear that we can perform a partial gauge fixing 
by requiring that (f>^ lies within the Cartan subalgebra and has the form (j^) of 
the matrix D above. This particular form corresponds to a fixed Weyl chamber 
within the Cartan subalgebra. The residual gauge group is then U{1)^^^. This 
gauge fixing procedure is described in details in Following we adopt the 
convention that the indices k, . . . are summed from 1 to — 1 (the indices 
of Cartan subalgebra) , while capital indices I,J,K,... from the middle of the 
alphabet run from N to N"^ — 1 (the indices of the subspace spanned by the 
roots). Also we assume that the indices a, 6, c, . . . run from 1 to 8 and correspond 
to the first eight coordinates in . Finally, let us adopt a new notation for 
the bosonic Cartan variables D'^ = t/)^ (that will be convenient in the future 
when we split the variables into the "fast" and "slow" ones). The gauge fixed 
Hamiltonian is H = Hq + Hb + Hp + H4 , where 
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Hn = -- 
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1 



Hp 



2 90ja(/)J ' 2 
i 

~ ^2 



(4) 

(5) 
(6) 



Hi 



fAIjfAKL4>'i<l^''AK^l 



-jAijlAkhDIDf^ 



f Aij f AK l{D°; (, 



1 
2 

'-flABri^AV'KB + ]^{w'w)ijLiLj 



(7) 
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The following notations arc used above: 

( d i 
Li = -fiBC y-'^B-Q^ + 2^BaAcc 

This Hamiltonian is self-adjoint with respect to the measure fl^ ^ d(j)\ Hi 
(the Faddeev-Popov Determinant was eliminated by redefining the Hilbert space). 

The terms Hb and Hp are Hamiltonians of bosonic and fermionic harmonic 
oscillators. To diagonalize these Hamiltonians, note that the eigenvectors of 
the matrices z'^j are the complex root vectors (m, n = 1, . . . ,N ,m ^ n). 
These eigenvectors satisfy 

z'ijEi^ = i(A^^. - \^!i)El^ (8) 

i^'mn)* = ^nm- (9) 

Here are eigenvalues of matrices D'^ = Y^iD'^T^. The vectors E^n define 
an orthonormal basis of a root subspace: 

Y,{El^J-El^ = 6^p6^, {m^n,p^q) (10) 
/ 

J2 «„)*<n = SlJ- (11) 

Note the equality Emn = E^^^T^ = -^emn , where emn is N x N matrix with 1 
on the m, n-th place and zeros elsewhere. These matrices satisfy the following 
commutation relations 

The eigenvalues of Hb and Hf (as we will see further) are 



Thus, we can introduce new variables <;i^„, A™", m,n = 1, ... ,N, m ^ n so 
that 

ri=J2'^^mnELn, A,<,= ^Ari?L- (14) 

The commutation relations for fermions are now 

{AJJ'",AP«} = (5c/35"f5'"« 
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and it is natural to set (A™")''^ — A"™ m < n. However, to diagonalize Hp 
some more work needs to be done. Namely, one should use the spm(9)-rotated 
fermions 

- A/3 rn<n 

\mn P"™A™"»n \ r, 

A„ = H^^ 7v^ m > n 

where 

jymn _ ^"in + {^m ~ ^Tt)^ ^ ^ „ 

"'^ ~ Yr Tr +A -A ) 

^1 mn\i mn ^ ^m ^n ) 

is an orthogonal matrix. Here and everywhere A„ with the suppressed upper 
index stands for A^. In terms of A™" variables Hp can be written as 

Hf=Y^ rrnn ((A™)1'A™" + A""(A"")t - s) (15) 

m < n 

where A± denotes the chiral components taken with respect to T^. 

Using (H) - we can rewrite all parts of the Hamiltonian in terms of the 
variables Df , Ai^, 0° , A™". The corresponding expressions for Hq and Hp 
are given in (jj) and (^5|) . The Hamiltonian for bosonic oscillators reads now as 

m<n ^ ii'-i^ \T rnn/ ni<n 

As one can easily see from (|l^) and ( |T6| ) the ground state energies for the bosonic 
and fermionic oscillators precisely cancel each other (we have 8 bosonic modes 
for each Vmn ,m < n). The expression for H4 in terms of the new variables is 
rather long. We relegate it to the appendix A in order not to complicate the 
main text. 

The normalized state vector of the oscillators ground state has the following 
form 

|o) = |Ob)|Of) 

|0b) = [ n ^™n2V-M expl-Y, r^nrmn^nnA 
\m<n / \ m<n / 

= n ri(A":)V«cfc> (17) 

m<n a — 1 

where \Fock) denotes the Fock vacuum for the fermions A'"". It is not hard to 
check that ( p^ is invariant under the residual U{1)^^^ gauge transformations. 
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3 Born-Oppenheimer approximation and pertur- 
bation theory 

The important feature of Matrix model that makes the existence of a threshold 
bound state possible is flat directions of the potential. Under the gauge-fixing 
condition that we employed, the coordinates along the flat directions are Cartan 
variables D^. We group them together with their fermionic counterparts A^. 
The coordinates in transverse directions are along with their superpartners 
A™". Following the terminology of Born-Oppenheimer approximation we call 
the variables Df, A, "slow" and the variables (?!)^„ , A™" "fast". When slow 
variables are considered to be fixed the dynamics of the fast ones is governed 
by the oscillator Hamiltonians Hb + Hp and by H^. In the asymptotic region 
where the frequencies r„i„ of the oscillators are large and assuming that 
can be treated as a perturbation, it is natural to expect that the fast degrees 
of freedom will remain in the oscillators ground state. The right assymptotic 
region for our purposes turns out to be Df oo . This is the same as setting 
Am — An — !■ cx) for all TO < n . We will assume that Am — A^ , to < n are all of 
the order r where r — > oo. Taking into account that the variables (/>J^„ in the 
oscillators ground state are of the order 1 / ^Vmn one can estimate that is 
of the order 1 / -^/r , that indeed allows one to treat it as a perturbation when 
r ~f oo. 

Thus, we search for an approximate solution to the spectral problem 

{H -E)\^') ^ 

in the form j^*) = \-)\'^{Di, Ai)) where |-) is the ground state of superoscilla- 
tors. The general formalism for this problem was developed in We refer 
the reader for the detailed explanation of such generalized Born-Oppenheimer 
approximation to that paper. Here we just want to explain the main idea of 
the formalism and develop a perturbative expansion suitable for the problem at 
hand. If one introduces a pair of projection operators P = |-)(-|, Q ~ 1 — -P, the 
Schrodinger equation breaks into a system of two equations 

P{H - E)P\^) + P{H - E)Q\^) = 

Q{H - E)P\-^) + Q{H - E)Q\-ii) = 0. (18) 

The second equation can be formally solved as 

gi*) = -{Q{H - E)Q)-'Q{H - E)P\^). (19) 

Substituting this solution into the first equation, we get 

[P{H - E)P - P[H - E)Q{Q{H - E)Q)-^Q{H - E)]P\^) , (20) 

i.e. the "reduced" Schrodinger equation for = \-)\^{Di,ki)). The first 

term in (pG) is the conventional effective Hamiltonian for the Born-Oppenheimer 
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approximation. The second one constitutes the correction term. Now note that 
in the problem at hand |-) is the ground state of superoscillators that has a zero 
energy. If we split the total Hamiltonian a,s H = Hose + H' where Hose is the 
superoscillators Hamiltonian and H' is all the rest, then HoscY) = 0. A direct 
analysis of formulae (0)-(0) shows that QHoscQ — Q{Hb + Hf)Q scales like 
r and QH'Q scales like 0(1) (because of the part of Hq that depends on 
variables). This is the reason to treat QH'Q as a perturbation. Now we can 
apply a perturbation theory in QH'Q to the term {Q{H — E)Q)^^ in (^ ) : 

[Q{H - E)Q)-^ = {QHoscQ + QH'Q - E)-^ = 

^ ^ -QH'Q ^^ I, ^ + .-. (21) 



QHoscQ — E QHoscQ — E QHoscQ — E 

Here Q{H — E)Q is understood as an operator on the Q-projected subspace 
of the whole Hilbert space and thus operator inverse makes sense (by the same 
reason we can write E instead of EQ in this equation) . Substituting this pertur- 
bation expansion into reduced Schrodinger equation (|2^), we get the following 
expression for the effective Hamiltonian 

+ {-\H'Q- \- QH'Q- \- QH'\-) + ... (22) 

QHoscQ — E QHoscQ — E 

The "propagator" ^q-e ^^^^l^s like 1 /r and the last expansion can be used 
to compute the effective Hamiltonian to any desired order in l/r. Although 
formally H' scales like 0(1) because of the Hq term, the contribution of the 
order comes solely from the first two terms in (|2^). This happens due to 
the fact that -ffo|-) scales like l/r (see appendix B). Hence, the correction term 
we need to have the effective Hamiltonian up to the order 1 /r^ is 

- i-lH.Q \ p Qg4|-). (23) 

^tlosc^ — ^ 

Another way to get the correction term is using an ansatz method similar 
to the one developed in [^. Both ( ^ ) and the ansatz method give the same 
result. 



4 Calculation 

First we need to compute the main contribution {■\H\-) = {-IHqI-) + {■\H4\-) 
up to the second order in l/r. For calculations it is convenient to express 
Hq in terms of variables. We have a linear correspondence between the 
variables ZJf = A'^X!^. Since by definition J^i^i"^^ — "i-J^i^i^a the condition 
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{D^',D'') = -2tr = ZJfLif gives A(^, = jAj^i'f. Therefore 

H ^ 

where the last operator is considered as a restriction to the subspace of functions 
of A^, . . . , A^ which are annihilated by (which means that they depend 

only on differences X!^ — A^ ). Both bosonic and fermionic oscillators ground 
state vectors depend on X!^ variables. This should be taken into account when 
calculating averages of differential operators such as Hq. Some useful formulae 
for computing such averages are given in appendix B. After this preliminary 
work the direct computation yields 

^ * m<n 

In Hi we have terms of the order 1/r and 1 / r'^/^ but it turns out that they 
contribute only in higher orders than . The remaining terms give the answer 

^ , , ' pm ' mn 

p^m rriytn ^ 



4 E 



2 

p^ra^n^p 



2(Ayj^ A^)(Ap ^m) ^pm ~l~ 



2 _L ^2 

mn 



'^mn ^pm ( Ap A-^, ) ^ '^mn '^pm ( Ap Xn ) ^ 

The correction term that gives a contribution of the order 1 /r^ is 



(25) 
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p^m^n^p 

m<n '^''^ 



(An — Ani)(Ap — An) 



''^mpij'pm H~ '^mn ^^ ^np) ^mp'^pn'^nmij'mp 



(26) 



Evidently the last term in (g6|) cancels with the analogous terms from ( 24 ) and 
(psi). Those terms are similar to the ones arising in SU{2) computation (see 
The remaining terms look, at the first glance, as if they can hardly cancel 
each other. To see that this indeed happens we rewrite them in terms of A„ 
variables only, using the fact that 

Tmn = A„, - A„ + 01^1 ,m<n. 



Then (^5|) contributes 
7 



E 

rn<p<n 



,{Xm ~ Xp){Xm ~ Xn) (Am — Ap)(Aj, — A,i) (Ap — A„)(Am — Xn) ^ 
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The same procedure carried on the first two terms in (|2^) yields 

1 1 



16 E 



m<p<n 



(Am — Ap)(Am — An) (Ap — An)(Am — A„) 



One readily checks that the sum of these two contributions vanishes. Hence, 
the outcome of our computation is given by the formula 

2 dD^dD^ ■ 

In appendix B we collected some formulae which we found useful for the 
computation described above. 



5 Discussion 

First we would like to explain how the present results can be compared with 
those obtained in |^ for the SU (2) gauge group. The authors of use a gauge 
invariant variable R to specify the asymptotic region. Namely, R is one of the 
eigenvalues of matrix 

Here lower a and h are SU{2) indices that run from 1 two 3. Without loss of 
generality we can assume that index 3 corresponds to the Cartan subalgebra 
and the basis is such that the matrix <I> is diagonal (any element of SU{2) 
can be taken as the Cartan generator). Then, clearly, 

* 
$=10* 



where — 4''^4>'i, coincides with r^2 in our old notations. Hence, in this basis 
R coincides with r. 

The result of |^ for the effective Hamiltonian (of radial degrees of freedom) 

is 

1 <f ^ d A 
"^"-~2dI^-RR-R^ 

while the measure on Hilbert space is R^^dR. Composing this operator with R 
from the left and with R~^ from the right we get the radial part of the Laplacian 
in 9 dimensions with the radial measure R^dR. Therefore, the radial parts of 
our effective Hamiltonian and the one found by Halpern and Schwartz are the 
same. It is not hard to check that the angular dependence is the same as well. 

Once we have the effective Hamiltonian, we can find the asymptotic ex- 
pression for Py^) = \-)\^{Di, Ki)) and then, using ( pj| ) and (pT[), one can get 
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the asymptotic form of the whole state vector \'^). Asymptotic sohitions to 
Heff\^{Di,Ai)) = have a basis of the form 

D-'-'^Yi,{D'^) ..... D-^'_-^"-'Yi,_,{D^^_,)\A,, . . . , A^_i) 

where Yj. are S0(9) spherical harmonics. Further constraints of the super- 
symmetry, S0(9)-invariance and invariance with respect to the Weyl group of 
SU{N) (i.e. the permutation group Sjy) should be put on A.;)) to single 

out the candidates for the ground state. 

Another comment we would like to make here is about the computation of 
Witten index. In paper ^ it was shown that Witten index for the model at 
hand has two contributions none of which can be made vanishing by the choice of 
the limiting procedure. The following expression for the correction (or surface) 
term was found 

Isurf = f treni-lfQW. 

^ JNf(R) 

The main ingredient in this formula is W - the approximate Green's function 
in the limit of large separation. The authors of ^ showed that one can take a 
free propagator for W in the SU{2) theory and gave the complete calculation 
of index in this case. In JT^ the surface term was calculated under a similar 
assumption for W in SU{N) theory. The present paper justifies this assumption 
and therefore completes the proof that the full Witten index for the SU (N) case 
is equal to 1. 
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A 

In this appendix we give expressions for the summands appearing in in terms 
of D'^ A (h°- A™" 

(Al) 

\fM.jfAKL{D'}4>''j - D'.rM^l = 
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hw*w)ijLiLj = — ^ .^ LnmLmn (^-5) 

^ — — Am, J 



where 



p 



B 



When working with bosonic oscillators it is convenient to introduce the creation 
and annihilation operators as follows 

V mn 



[«Ln, (flp,)^] = S'^'S^j^Sn, m ^ n,p ^ q. {B.l) 

Since both bosonic and fermionic oscillators ground states depend on eigenvalues 

, a certain care needs to be taken when calculating the averages of differential 
operators over the oscillators. The following formulae are useful 



{ s:mk ^nk\ 
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Fermionic quadratic correlator is 



1 '^^^(A'"")tr''A™"]|-). (BA) 



^ V / a/3 

for any m ^ n. 

Note that spm(9)-rotated fermions A™" depend on AJJj. The following cor- 
relator thus appears in computations 



(OF|^(A™")tr°A™"|Oj.) = --rr(r9[g"" ^^^^ ^ ,n). (i?.6) 
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